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ABSTRACT 

The flow in a two-dimensional curved channel driven by an azimuthal 
pressure gradient can become linearly unstable due to axisymmetric perturba- 
tions and/or nonaxisymmetric perturbations depending on the curvature of the 
channel and the Reynolds number. For a particular small value of curvature, 
the critical Reynolds number for both these perturbations becomes identical. 
In the neighborhood of this curvature value and critical Reynolds number, non- 
linear interactions occur between these perturbations. The Stuart-Watson 
approach is used to derive two coupled Landau equations for the amplitudes of 
these perturbations. The stability of the various possible states of these 
perturbations is shown through bifurcation diagrams. Emphasis is given to 
those cases which have relevance to external flows. 
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1 . INTRODUCTION 


In plane channel flow, instability arises due to the amplification of 
Tollmien-Schlichting (TS) waves. As these waves grow, they modify the mean 
flow, produce higher harmonics, interact with other waves, and probably pro- 
duce turbulence. The initial stage of development of these waves from the 
linear region to the weakly nonlinear domain was analyzed by Stuart (1958) who 
derived the Landau equation for the temporal development of the TS wave. The 
presence of the cubic nonlinearity in this equation modifies the otherwise 
exponential variation inherent in a linear theory. This theory was able to 
explain the existence of an equilibrium finite amplitude perturbation in cer- 
tain regions near the neutral curve. 

Taylor (1923) was the first to consider the instabilities that arise due 
to the curvature of streamlines. He investigated the flow between two concen- 
tric cylinders due to the rotation of the inner cylinder with the outer cylin- 
der stationary. He found that the flow becomes unstable when the parameter 
Re(d/Rj)l/^ (now referred to as the Taylor number) exceeds a value of about 
41. Here R^ is the radius of the inner cylinder, d (<<Rj) is the gap width 
of the cylinders, and Re is the Reynolds number based on the speed of the 
inner cylinder and d. The instability that appears as the speed of the inner 
cylinder exceeds the critical value is in the form of toroidal vortices. 
These vortices are modelled theoretically by an axisymmetric perturbation and 
they are stationary when they first appear. 

Dean (1928) also investigated the instability in a curved channel due to 
the curved streamlines (see Figure 1). The flow in his experiment was gener- 
ated by an azimuthal pressure gradient. The channel is formed by portions of 
two concentric cylinders having channel width d « R^ . Basing the Reynolds 
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number, Re, on the mean speed of the unperturbed flow, Dean (1928) and 
Walowit, Tsao, and DiPrima (1964) found that instability arises when 
Re(d/R^)*^ exceeds a value of about 36* Here, too, as in Taylor's experi- 
ment only axisymmetric disturbances were considered. 

In a detailed analysis of the linear stability of curved channel flow, 

Gibson and Cook (1974) argued that in a curved channel of very small curvature 

non-axisymmetric disturbances can play a significant role in destabilizing the 

mean flow. Such perturbations are analogous to TS waves in a plane channel. 

Their linear stability analysis shows that for channels with very small 

curvature, the critical Reynolds number for the TS waves is almost independent 

of n , (n « R^/R^, R^ = radius of outer wall), and it approximates very 

closely the corresponding value for a plane channel. The critical Reynolds 

number for the axisymmetric instability (Gortler vortices), on the other hand, 

is quite sensitive to n for n close to 1 (Figure 2). For a particular 

value of n = n , the critical Reynolds numbers for these instabilities are 
c 

identical. For a slightly wider channel, the critical Reynolds number for the 
Gortler instability is lower than the almost constant critical Reynolds number 
for the TS perturbation. For a narrower channel, the critical Reynolds number 
for the Gortler instability is higher. It is, therefore, reasonable to expect 
that near both perturbations could exist simultaneously and thereby 

interact with each other. 

The purpose of this paper is to analyze the weakly nonlinear interaction 
of these two instabilities (one axisymmetric and the other non-axisymmetric) 
which arise when the radius ratio is nearly We use a multiple scale 

version of the Stuart-Watson method approach to derive the two coupled 
ordinary differential equations for the amplitudes of these perturbations. 
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While these two equations cannot be solved explicitly, they nevertheless yield 
significant information about the various possible bifurcations that can take 
place in the presence of these perturbations. Moreover, the stability prop- 
erties of the equilibrium states can be deduced. 

When considering a growing boundary layer, a self-consistent analysis of 
wave interactions within it requires the application of the triple deck theory 
as shown by Hall and Smith (1984). In channel flows, however, we do not need 
to consider the effects of boundary layer growth, thereby greatly simplifying 
the analysis while still giving a qualitative picture of what might happen in 
an unbounded flow. It is in this context that we wish to study the Gortler/TS 
interaction in a curved channel. This study may be viewed as an extension of 
the work of Gibson and Cook (1974) into the weakly nonlinear regime. 


2. MEAN FLOW AND PERTURBATION EQUATIONS 

Let (r, 0 , z) be the cylinderical coordinates with the axis of the 
concentric walls along the z-axis and and R Q the radii of the inner and 
outer cylinder respectively (see Figure 1). 

When the flow between the concentric walls is maintained by a constant 
azimuthal pressure gradient yg- (< 0), the solution of the momentum equa- 
tions yields 

U(r) = W(r) - 0, (2.1) 


where U(r) and W(r) are the radial and axial mean velocities, respective- 
ly. The azimuthal velocity is given explicitly by 
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R 2 

„/ x O ap r . n log n ,Jl n i 

V(r) = _ _ r log r + 5 - (r - 1)J , 

2vp 30 r( 1 -ri ) 


( 2 . 2 ) 


(n < r < 1 ) 


where v and p are the kinematic viscosity and density, respectively, 

and n = R /R • Here the distance from the axis is normalized with respect 
i o 

to R q , so that r varies from n at the inner radius to 1 at the outer 
radius • 

A channel with small curvature will behave locally like a plane channel, 
and for it the azimuthal velocity V(r) should approach the familiar para- 
bolic shape. As n approaches 1, the velocity 


v(r> ■ - 2 ^!f [° - 11)2 1(1 - 5> ! (2 - 3) 

= v c(l - c), 0 < c < 1 (2.4) 

m — — 

£ 

where V m = (1 " n)2 (2 * 5) 

and £ is given by 

r - (1 - n)S + n • (2.6) 


Thus ? varies from 0 to 1 as r varies from n to 1. Note that 

in this limit, V m is four times the center line velocity. 

Based on (2.3), (2.4), (2.5), and (2.6), (2.2) can be simply expressed as 

V(r) = V f(r) 
m 


(2.7) 


where 


f(r) = - — - — y [r log r + — g ? n (r 2 - 1)] 

(l-n) r(l-n ) 


( 2 . 8 ) 


The fully nonlinear disturbance equations for the radial velocity u, 
azimuthal velocity v, axial velocity w, and the continuity equation are as 
follows: 


5? " " W ~ F VST ' 7T "z'^e 


3 p (l-n) 3 ,3v. 1 /3wv 

~ " '^2 

(!-n)9 


, r t (l-n) . 1 ~ f(r) (l-n) . \ 

+ l Ri + 0 9 + — 9 „ " T— 8 J U 


„ 2 2 90 _ 2 zz r Re 9- 

Re r Re 


+ (-■^2!^ S„+4 £(D s e b 


Re 2 9 r 

r 


(2.9) 


-a-n)[^ TV -^ + 7T77-l"lf + “V)] • 

Re (l-n) 


8 v Re v. 1 9v 

7 r ~ T~ (1 " n)a e p ~ 7 a “ ^ ac 


♦[-fe 


df ( r) 


f(r)- 


+ [te(l - I)9 t + 9 99 - > M + Re iiii (1 - n)3 


f (r) 


( 2 . 10 ) 


+ Kr (1 - n) 2 ] 


+ (1 ’ n) [ir^7ff + Re 7 V + r + u^yVl • 


uv . w 


6 - 


2 

9 w „ 2 „ 
— 2 = Re 3 z p 

3 ? 


( l -n ) 9 w 
r 9£ 


[Re(l - 


(1-n) „ 

n)9 t 2 8 ee 

r 


8 zz ♦ Red - TO ^3„]» 


( 2 . 11 ) 


+ (1 - n)Re — 9 w + u + w9 w , 
r o d c, z 


9u _ (1-n) Re /, 

9? r u “7" (1 " n)3 0 v " 9 Z W 


( 2 . 12 ) 


where the quantities that appear in (2.9), (2.10), (2.11) and (2.12) are non- 
dimensionalized versions of primed physical quantities shown below: 


Pressure d = — 


Radial position r = r'/R ** (1 - n)C + n> n .< r ^ 1 


0 < c < 1 


Axial position z = z'/d 


Azimuthal position 0 = 0' 


Time t - t'/(=2) 


Radial velocity u » u'/(—) 


Axial velocity w = w'/(t0 

d 
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Azimuthal velocity v = vVV m * 

Note that u^ and w' are scaled with respect to the diffusive velocity 
scale while v" is scaled with respect to the convective velocity scale. The 
mean flow is (0, V m f(r), 0) and the nondimensionalized perturbation is (u, 
v, w) in the (r, 0, z) direction. The Reynolds number 

Re - V d/v, 
m 9 

where d = R Q - Rj_. 

The effect of the purely azimuthal (non-axisymmetric) and the purely 

axial (axisymmetric) perturbation can be modelled by a general expression for 

the perturbation proportional to exp [at + i(kz + m0)] where the non- 

dimensional axial wave number k = k'd, the non-dimensional azimuthal wave 

number m - m% and the non-dimensional complex growth rate 
V 

a = or'/(-^) = a + i a . The partial derivatives a 9 8 f 8 can now be 
R R X t z 0 

O 

replaced by a, ik, im respectively within the linear part of the equations 
(2.9) to (2.12). 

The four equations of motion can be written as a set of six first-order 
ordinary differential equations, as was done by Eagles (1971). Thus we write 

(P, , u » V, w) T = q (2.13) 

so that the equations can be written as 

|jr q_ = * q. + n 


(2.14) 
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where * is a 6x6 matrix representing the linear contribution and jn 
is a six element vector containing the nonlinear terms. The last three ele- 
ments of the vectors q_ and JL are zero at t ^ ie two channel walls. The terms 
representing the linear and nonlinear contributions will be discussed in the 
next section after an explicit expression for q is given. 


3. PERTURBATION EXPANSION FOR NONLINEAR WAVE INTERACTION 

The perturbation to the mean flow is expressed In the form: 

q = e{(A E + B F) + c.c.} 


+ e 2 {(C E 2 + D F 2 + G EF*+ H EF) + c.c.} 


+ e 2 { J E° + K F°} 


3 3 3 2 * 2 

+ e { (L E+MF+NEF +PEF 


+ Q EF 2 + R E*F 2 ) + c.c.} 


+ e 3 {(S E + T F) + c.c.} + 0(e 4 ), 


(3.1) 


where _q_ is the perturbation vector given by (2,13), c.c, represents the com- 
plex conjugate of the terms in the preceding parentheses, e is a small ex- 


pansion parameter, 


represents complex conjugation 
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The vector 


r, 

A = ( a i » a 2» a 3» a 4» a 5> a 6^ 


represents the Gortler perturbation, 


A - (*>1 > b 2 » l 5 3 > l 5 4 » b 5 » b 6 ^ T 


represents the TS perturbation, 


(3.2) 


(3.3) 


and 


E = exp(ikz), 


F = exp(im0 )exp(ot) , 


(3.4a,b) 


a = a R + icj . 


(3.5) 


A E and B_ F represent the Gortler and TS instabilities respectively. 
Other terms represent higher harmonics and the mean flow modification which 
arise due to the quadratic nonlinearity in the equations. 

The amplitude coefficients £, D_, Gj, and _H are due to the direct 
interaction between the TS and the Gortler instabilities. As self interaction 
occurs for each of these instabilities, the mean flow profile itself becomes 
modified through the generation of effects represented by J_ and K_. The 

o 2 

terms at order e arise due to the interactions of the terras of order e 
•• 

and the Gortler and TS perturbations. 

As will become clear later, it will not be necessary to solve for the 

O 

unknown coefficients at order e • By using the solvability condition for 
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the equations governing j>_ and T_, the evolution equations for the amplitude 

M 

of the Gortler (A_) and of the TS (B_) waves can be found. The method of 
multiple-scales used in this paper for obtaining these evolution equations 
follows closely that suggested by Matkowsky (1970). 

Throughout, we will require that the last three components of the vec- 
tors A to T are zero at the two curved walls in order to satisfy the zero 
velocity boundary conditions. 

For the perturbation expansion (3.1), the appropriate slower time vari- 
able is t defined by 

t = e 2 t, (3.6) 


where e is the same small expansion parameter used earlier. This gives 


3 

3 1 




3 

3t 


+ e 


2 3 
3t * 


(3.7) 


The Reynolds number Re 
bility curve so that 

Hence, 

and 


is expanded about a value 
Re = R Q + e 2 R x . 

Rl 2 

1/Re « 1/R n - -j e 
R 0 

~ 2 ? 
1/Re « 1 /Rq - 2 — j e . 

R 0 


Rq on the neutral sta- 

(3.8) 

(3.9) 

( 3 . 10 ) 


Returning to (3.1), we now regard all the vectors to as functions 
of both £ and t. For example: 
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A -► A U) X(x ) 

B «• B(C) Y(t ) 

which states that in the neighborhood of the neutral stability curve the solu- 
tions for the linear problem, A(5) and B(c)> adequately represent the 
shape of the perturbations in the weakly nonlinear regime. The amplitudes of 
these perturbations, X(x ) and Y(t), depend on the slowly varying time 
variable x . 

The remaining vectors can be represented similarly, e.g., 

C > C(S) X 2 (x) 

H H(c) X(x) Y(x). 

Henceforth, when referring to vectors A. to JC we will only be con- 
sidering their spatial dependence because it can be shown that the temporal 
dependence cancels out throughout the equations for these vectors. 

By substituting (3.1) and (3.7) to (3.10) into (2.13), we obtain an equa- 
tion for q_ of the form: 

+ £ 2 L 2 [9 z ,3 9 ] a 

+ e 2T 2i[B t ] i + e 2T 22 t 8 x ] 3. (3.11) 

+ N 1 [3 z ,3 0 ] {£® q} 

+ e 2 N 2 [3 z , 3 0 ] {q® q} . 
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1 * 1 , L 2 > Tp T 2 i» T 22 are linear operators, equivalent to 6 x 6 matrices. 

Their dependence on the differential operators 3^, 3^ , 3 fc and 3 is 

indicated within the square brackets. Both ^ ^ z * ^ 0 ^ ^3.® 3) and 

N„[3 ,3 fl ] {q®q} represent nonlinear terms containing the operators 3 

and 3^ . (q (x) 5 } symbolizes quadratic terms comprised of the components 

of Here onwards we shall refer to — l^z’^9^ ^3.® 3.} by . Since 

3 

in the analysis of (3*11) we will only be concerned with terms of 0(e), it 

2 

will not be necessary to consider the term e (q® q} which is 

4 

0(e )• The elements of the operators are given explicitly in Appendix A. 

•• 

To determine the linear stability problem for the Gortler and the TS per- 
turbations acting individually, terms of 0(e) need to be considered. From 

these terms, those with coefficient E will give the linear stability pertur- 

• • 

bation equations for the Gortler disturbance. Similarly, the terms with 
coefficient F will give the equations for the TS disturbance. 

At 0(e) 

(i) The terms with coefficient E give 

dA 

3 - - I^fik'OlA + TjfOlA (3.12) 

for the Gortler perturbation. 

(ii) The terms with coefficient F give 

dB 

■%r = LjlO.imlB + Tjlio^B 


(3.13) 
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for the Tollmien-Schlichting perturbation. 

From these equations, vectors A_ and B_ can be determined. Note the 

difference in the parameters of the operators in (i) and (ii), particularly 

the fact that 3 fc ■ 0 for the Gortler while for the TS, 3 * ia^. In 

both cases we are considering neutrally stable perturbations. 

2 

When collecting terms of 0(e ) for determining vectors _C to K_, the 

operators 1 * 2 , T 21 , T 22 can be neglected because when they operate on _q_, 

2 

contributions of 0(e ) are produced. As these operators are already pre- 
multiplied by e^, the net contribution from these terms will be O(e^) 
and hence negligible. In the equations for C_ to K_ that follow, the non- 
linear contributions are contained in _Nj. Of the six elements of the vector 
representing the contribution from Nj , only the first three (nj, n 2 , ng) are 
non-zero, and these are listed following each equation. The discussion 
following (3.5) gives the physical basis for presence of the vectors C_ to 
K. 


At 0(e 2 ) 


(i) The terms with coefficient give 


dC 

= L 1 [2ik,0]£+ Tj [0]£ + nonlinear contribution from , (3.14) 


n, = — (1 - T|){ — — 5 a c/ r } 

R o r 


where 
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n 2 “ a 4 a 2 + a 4 a 5^ ” + ik a 6 a 5 

2 

n 0 = a.a 0 + ik a, . 

3 4 3 6 


(ii) The terms with coefficient give 


d D 

— = L^[0,2im]]} + T^[2 10^.]]) + nonlinear contribution from , (3.15) 


where 


n = (1 - 


n 0 = (1 - 


2 2 
bf bf 

"> tr- r* 


V 

2 

b.b 0 imb c R ri b.b. 


= (1 - n) — imb^bg + b^b^. 


(iii) The terms with coefficient EF* give 
dG 

= L^[ik,-im]G + Tj[-iOj]£ + nonlinear contribution from N ^ , (3.16) 


where 


2a. b * 

n x = (1 - n) { 2 

R 0 r 


* 

im a,b c ik b. a, 

4 5 4 6 


R O r Rq(I-ti) Rq(1-tO 


ii * ii 

ik b. a , im a c b, 2a c b c 

6 4 5 4 5 5 i 

r J 


R o r 


n 2 = (1 - ti) { 


a 4 b 2 + b 4 a 2 


( l-n) 


im a c b_ a.b.+b.a- ik a_b. 
055^4545^ 5 6 1 

r r (l-n) 1 


3 z, v * * * * 

n 3 = -im(l - n) R q a 5 b 6 + a^ + b^ + ik b^ . 
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(iv) The terms with coefficient EF give 
dH 

= Lj[ik,im]Il + Tj[io^.]Il + nonlinear contribution from , (3.17) 

where 


n, = -(1 - n) 


—S— 2a, b, - a. im b — - b.ika, + — s— ^ b,ika 

R 2 r 4 4 l^r 4 5 R 2 ( ,_ n) 4 6 R 2 (1 _ n) 6 4 


0 

a c im b. 
5 4 

V 


O' 

2 a 5 b 5 


= 


„ „ ,‘4 b 2 * b 4 a 2 . V” b 5 t a 4 b 5 + \ a 5 . V ka 5 , 

R o (1 • 1 yn) + — f — + — iy + (l-n)R 0 t 


n„ = 


R 0 (l-n) 

im a c b, + a. b. + b. a_ + ik a,b, . 

r 56 43 43 66 


(v) The terms with coefficient E® give 


dJ 

d? 


= Lj[0,0]J[+ Tj [0] J[ + nonlinear contribution from , (3.18) 


where 


2(l-n) * . 2ik f * *1 . 2(l-n) 

= — a. a. + — = — ]a.a, — a.a, f + 

R»r 4 4 „2 1 4 6 4 6 J r 


0 


a 5 a 5 


0 


* * 


* * 


(l-r\ ) * * 

a 4 a 2 + a 4 a 2 + — r — (a 4 a 5 + a 4 a 5> + lk(a 5 a 6 “ a 5 a 6 ) 


n 3 = a 4 a 3 + a 4 a 3 * 


(vi) The terms with coefficient give 
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dK 

dif 


L^OjOlK + [ 0 ]K + nonlinear contribution from Nj , (3.19) 


where 


b 4 b 4 . im 


n l “ 2(1 ^~2 + R~F (b 4 b 5 ‘ b 4V + b 5 b 5^ 


R o r 


(l-n) 


n 2 = b 4 b 2 + b 4 b 2 + f b 4 b 5 + b 4 b 5 > 


n 3 


R 0 (l-n) 


r { imb 5 b 6 “ imb 5 b 6^ + b 4 b 3 + b 4 b 3 * 


So far we have collected terms with coefficient e and obtained the 

equations for the spatial dependence of the neutrally stable Gortler and TS 

perturbations followed by collecting terms that have coefficient e^. At 
3 

0(e ), we will see that only equations of _S_ and T_ need to be considered 
for obtaining the time evolution of the Gortler and TS amplitudes. S_ and 
X are functions of c and x . In the following, we shall again revert to 
writing the Gortler and TS perturbations as A X(x) and Ji Y(x) respec- 
tively, with A_ and B_ dependent only on 5 . It will be seen from the 
form of the equations for S_ and T_ that the temporal dependence does not 
cancel out; indeed, it is this very property that allows us to get time evolu- 
tion equations for X(x) and Y(x). 

O 

Collecting terms of 0(e ) with coefficient E, and with coefficient F, 
we obtain the equations for S and T: 

9_S 

jj-T = L^ik.O]^ + T^O]^ 


+ L 2 [ik,0]A X(x ) + T 21 [0]A X(t) + T^fa^jA X(x) 


(3.20) 
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and 


+ contribution from 
8T 

g^- = L^O.imjT + TjIiOjlT 

+ L 2 [0,im]B Y(t ) + T^riOjjB Y(x) + T^fa^jB Y(x) (3.21) 

+ contribution from . 


The contributions from to both these equations involve a very 

large number of terms and are therefore not written explicitly at this 
stage. They will however appear in the final equations for X(x) and Y(x). 

The homogeneous parts of (3.20) and (3.21) are the same as those of 
(3.12) and (3.13) for ^ and ^ respectively. In order that the non- 

homogeneous equations (3.20) and (3.21) have solutions, the nonhomogeneous 
parts of the these equations should be orthogonal to the adjoint column 
vectors and JB respectively. These vectors are solutions of the 

following equations: 


T 

= -[Ljfik.O] + TjfO]] 1 A (3.22) 

dB T 

•fir = -[LjfO.im] + T^iOj]] 1 B (3.23) 


^ rj% ^ ^ ^ ^ ^ ^ rji 

where _A = (a^ ,a 2 ,a 3 ,a^ ,a 5 ,a^) and IJ = (bj ,b 2 ,b^ ,b^ ,b 5 ,b^) . Unlike 
the boundary conditions for the equations for A_ and B_ where the last three 
of their components are zero, here the first three components of A and 


are zero at the boundaries 
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Using the orthogonality condition, we obtain the following equations 
for X(t) and Y(x): 

/ dc[A T L,[ik,0]A X(t) + A T • T [9 ]A X(t) 

0 

(3.24) 

+ A^ • (contribution from N^)] =0 

and 

/ d?[B T L [0,im]B Y(x ) + B T T [io ]B Y(t ) 

0 

(3.25) 

+ B^ T 22 ^t^— + — * (contribution from N^)] = 0 . 

The terms representing (A^ • contribution from ) and (B^ • contribu- 
tion from N^) are given in Appendix B. After integrating over 5 , these 
equations can be written as 

d §~ = R x Bj X(t ) + X(t)|x(t )| 2 + n x X(t)|y(t )| 2 (3.26) 

and 

d = Rj B 2 Y(t ) + 6 2 Y(x ) | X( t ) | 2 + n 2 Y(t ) I Y(t ) I 2 (3.27) 

where 3 ^, 6 ^, an ^ U = 1>2), are coefficients obtained from (3.24) and 

(3.25). Rj is a measure of the deviation from the neutral stability curve as 
is given by (3.8)* It appears in (3.26) and (3.27) because it is a common 
factor in matrix L 2 in (3.24) and (3.25). Equations (3.26) and (3.27) are 
the coupled Landau equations which determine the time evolution of the ampli- 
tudes of the Gortler and TS perturbations. The analysis of these equations 
will be presented after the next section. The following section gives a brief 
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description of the numerical method used for obtaining the coefficients _A 
to K_ in (3.1). All the terms in the perturbation expansion have been 
verified by using the symbolic manipulation language MACSYMA. 

4. COMPUTATION OF THE COEFFICIENTS 

The equations governing the 5 dependence of the coefficients _A to 
K are given in the previous section. A^ and _B_ are described by a set of 
homogeneous ordinary differential equations while the equations for the re- 
maining amplitudes are nonhomogeneous. 

A fourth-order finite difference scheme (Malik, Chuang, and Hussaini 
(1982)) was used to solve these equations. For details of the method, the 
reader may refer to Malik, et al. (1982) and Hall and Malik (1986). The cal- 
culations were performed on a nonuniform grid which clusters the points near 
the walls. A suitable distribution of grid points was obtained using the 
relation 

5^ = (sin(irx i /2) + l)/2, (4.1) 

where 

= -J-j- [2i - 1 - N] , 1 < i < N (4.2) 

and 

N = total number of grid points. 

In order to determine the vectors A^ to K_ to 3 digit accuracy, 51 grid 
points were sufficient for the range of Reynolds number and wave numbers that 
we considered. 
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The coefficients of the Landau equations derived in the last section are 

functions of the vectors _A to K. These vectors become dependent on the 

slow time variable t only when we perturb the flow from its neutrally 
stable state* In the neighborhood of the neutral stability curves for the 
Gortler and TS perturbations, the shape of the vectors (as functions of c) 
remains unchanged, and hence we only need to consider the shapes of the 

neutrally stable modes. 

The familiar neutral stability curves for the linear Gortler and TS 

perturbations are presented in Figures 3 and 4, respectively. Here, the 
generally accepted convention of labelling one of the arms of the stability 
curve as "lower" and the other as "upper" is used. To analyze the different 
kinds of possible interactions between a Gortler and a TS perturbation at an 
arbitrary Reynolds number Rq consider the schematic diagram of Figure 5. 
(Refer to this figure and its caption for the abbreviations GL, GU, TSL, TSU 
used in what follows.) A Gortler perturbation GL with wave number k^ and 
Reynolds number Re slightly different from Rq (Re = Rq + e R^ ) can 

interact with a TS wave with the same Reynolds number but with wave numbers 
mj or m 2 corresponding to TSL and TSU respectively. Similarly, GU with 

wave number k 2 can interact with either TSL or TSU. So, in all there are 
four possible interactions. 

It can be seen from Figure 2 that for n = R./R 0 = r\ = 2.179 x 10 

12 c 

the critical Reynolds number for the Gortler and TS perturbations is 
8 x 5772.2 - 46176 where 5772.2 is the critical Reynolds number for a 
plane channel flow based on half channel width and centerline velocity (Orszag 
(1971)). It is for this value of r\ that we compute the amplitudes A_ 

to K_ for 46176 _< R^ _< 120,000. This range is probably sufficient to 
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reveal the possible interactions. Results for other values of p in the 
neighborhood of r\ ^ can be obtained using a simple argument that we will 
present in the next section. 

The components of A to K for each value of Rq are used to compute 
the coefficients 0p 6j, n^, 3 2 > ^ 2 * atlc * ^2 t * ie tWO ^ anc * au equations. 
For computing these coefficients, vectors A_ and need to be normalized. 
This was done by dividing _A by its centerline azimuthal component (center- 
line value of and dividing _B_ by its centerline radial velocity com- 
ponent (centerline value of b^). Since 0^ and 6^ do not depend on the 
presence of a TS perturbation, each of these has a unique value for each point 
on the neutral stability curve. When both GL and GU are considered for a 
fixed Reynolds number, 0^ and 6^ will each have different values on the 
two arras of the neutral stability wave, corresponding to the different wave 
numbers. A similar argument applies to 8^ an< * ^2 are independent 
of the Gortler perturbation. These coefficients will be further discussed in 
the next section. 


5. SOLUTION OF THE LANDAU EQUATIONS 

In this section, we analyze the possible interactions by studying the 
properties of the coupled Landau equations. These properties are displayed in 
the form of bifurcation diagrams which show the amplitudes of the equilibrium 
states and their stability properties. 

Equations (3.26) and (3.27) can be written in terras of | X | 2 and | Y | 2 : 


. „ 6| „ |x | 2 + 6,„|x| 2 |x| 2 ♦ nJx| 2 |Y| 2 


1R 


1R' 


' 1R 1 


(5.1) 




2i„i2 


2 1 „ 1 2 


-U g 2R |Y|‘ + «,. |YHxr +n 0B |Y|‘|Y 


2R 


2R 


(5.2) 
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where S 1R , 3 2R , <S 1R , 6 2R , P 1R , and p 2R are real parts of the corresponding 
Landau coefficient, e.g., 3 2 = 3 2R + i8 2I » It is found that ^ is real 

and so 3^ = 3^ R . By suitably scaling the amplitudes X and Y it is 
possible to make 6^ = -1 and | n ^ I “ 1 so as to facilitate the analysis 

of the equations. For notational convenience, we replace Rj in (3.26) and 
(3.27) by y in the above equations. From here on we shall drop the sub- 
script R because all the coefficients of the equations are real. 

The growth rates with repect to t (for y = 1) of the Gortler (3^) 
and TS ($ 2 ) perturbations are shown in Figure 6 and Figure 7 respec- 
tively. In Figure 7 it should be noted that the negative values of 3 2 
corresponds to TSU in Figure 4. 

There are two graphs each for p^ and 6 2 depending on the types of 

possible interaction between a Gortler and TS wave. Figure 8a shows p^ 

versus Reynolds number for the interaction of GL with TSL and TSU and Figure 

8b displays the same variables for the interaction of GU with TSL and TSU. 

Figures 9a and 9b give graphs for 6 2 for the same interactions as given in 

Figures 8a and 8b. Figure 10 shows a graph of Reynolds number versus p 2 » 

The graphs of the coefficient of the Landau equation mentioned in the 

last two paragraphs have been computed for (1 - p) = (1 - p ) = 2.179 x 

10 \ which corresponds to a channel with very small curvature and one for 

which the critical Reynolds number for the Gortler and TS perturbations are 

identical. In what follows we will extend the analysis to channels for p 

in the neighborhood of p . 

c 

In (5.1) and (5.2) the coefficients of |xp and | Y | ^ are y3^ 

and y3 2 respectively, which shows that the linear growth rates of |xp 

and | Y | 2 are proportional to the deviation y from Rg (Reynolds number 
2 

= Rq + e y). Similarly, it is reasonable to assume a linear dependence of 
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the growth rate on the deviation from the radius ratio r^. If we write 

o 

(1 - n) = (1 - n ) + e v as a perturbation from (1 - n c )» and take e 

to be the expansion parameter given by (3.8), the growth rates of |x|^ and 

| Y | 2 would also be linearly dependent on v. The effect of the deviation 

from (1 - T) ) on the coefficients of the nonlinear terms of the Landau 
c 

equations is of a higher order than we are concerned with and so we shall only 
consider its effect on the growth rate through the parameter v. 

Equations (5.1) and (5.2) can now be written as 

y = Mi + v V ,x|2 + 6 i l x l 2 l x l 2 + n i l x l 2 M 2 (5 * 3) 

(y6 2 + v V l Y l 2 + 6 2 I y I 2 I x I 2 + n 2 I y I 2 I y I 2 * (5 * 4) 

The zero-growth rate curves for |Xp and |y| 2 are straight lines in 
the (v,u) plane, passing through the origin and having slopes of (-y^/f^) 
and (— y 2 2 ^ res P ect ^ ve ^y* Numerical computations show the first slope to 

be of order —10^ and the second to be approximately zero; therefore in what 
follows, we take ~ 0* To compute these slopes, the wave numbers for the 

Gortler and TS waves were fixed at their values for p * 0, v = 0, which is 
equivalent to R = Re c and (l-n)*(l-n c ). Th e slopes were found for 

this neighborhood for these fixed wave numbers. As can be expected, the 

slopes are of the same order of magnitude as those of the curves for the 
critical Reynolds number versus (1 - n) given by Gibson and Cook (see 
Figure 2). Note that the wave numbers change along their curves, while in our 
case we keep them constant. We computed the change in Reynolds number with 
change in (1 - n) at the cross-over point shown in Figure 2. 
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There are four possible steady state solutions to (5.3) and (5.4): 


(i) |X | 2 = 0, |Y| 2 = 0 

9 9 3 1 VY i 

(ii) | y | z = o, |x | 2 = - (p + ^ 

6 i e i 

(iii) | X | 2 = 0, | Y I 2 = - U 

n 2 

9 VY 1 

(iv) | x | = [-n 1 e 2 y ■*" n 2 6 j (y + -g — ] / Cn x 6 2 - 6^) 

9 vy i 

I Y I = [-^2^1^ + g — ^ + 6 i^2 y -l /^ n i <S 2 * <S 1 T1 2^ * 


(5.5a,b) 


(5.6a,b) 


(5.7a,b) 


(5.8a) 


(5.8b) 


To assess the stability of these states, it is necessary to linearize the 
equations about these states. The reader may refer to Boyce and DiPrima 
(1977) for a discussion of the stability analysis of such coupled equations. 

Table 1 summarizes the values of aTU * 5^,1 = 1,2 , for the 
four possible interactions. This tabulation is mainly to facilitate the 
analysis and discussion of these interactions; other values can be obtained 
from the graphs for these quantities. This table of results can be used to 
show that many possible equilibrium states exist depending on the Reynolds 
number. Here we shall concentrate on the three cases which we believe to be 
of most practical importance. The bifurcation pictures for the other cases 
can be found in say Keener (1976) or Guckenheimer and Holmes (1984). 

The three cases we consider are: 


(a) Interaction of Gortler and TS waves for 


n - 


n + e 2 v 

C 


and 
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2 

Re = Re c + e p . Here we consider the cases when the Reynolds number 
is at or very close to the critical value for both perturbations, and 
with wave numbers corresponding to the critical Reynolds number and 
its vicinity. Representative values for this case can be found near 
the end of each of Sections A - D in Table 1. 

(b) Interaction of TSL with GL. 

(c) Interaction of TSL with GU. 

Case (a): In this case, it is found that r^/r^ > ^ ^1^2’ 

6^ < 0, > 0, > 0 for i = 1,2. For v > 0, the Gortler mode is the 

most unstable on the basis of linear theory and for v < 0 the TS wave is 
the most unstable. 

The solution (5.5a,b) exists for all values of p whilst (5.6a,b) and 
(5.7a,b) exist for p > and p < 0 respectively. The mixed mode 

solution (5.8a,b) can exist for either a finite, zero or semi-infinite range 
of values of p depending on In the present case, we find 

that if v < 0, the mixed mode does not exist. However, for v > 0 the 
mixed mode exists for a finite range of values of p including the 
origin. The bifurcation diagrams for this case are shown in Figures 

(lla,b). In these figures, continuous and broken lines correspond to stable 
and unstable solutions of the Landau equations respectively. We note that the 
TS mode can never be in stable equilibrium without the presence of a Gortler 
mode. In contrast to this situation, the Gortler mode can exist alone and be 
stable to small perturbations. However, the finite amplitude states in these 
figures are unstable to sufficiently large perturbations. This instability 
leads to |x| and |y| terminating in a finite time singularity as in the 
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case for a TS wave in a straight channel. Thus, the threshold amplitude 
phenomenon of Meksyn and Stuart (1951) persists in the presence of Gortler 
vortices. It is not possible to quantify the effect of the Gortler mode on 
the threshold amplitude. However, a phase plane analysis of the Landau equa- 
tions shows that the Gortler mode significantly reduces the size of the finite 
amplitude perturbation required to induce the finite time breakdown of the 
equations. In that sense, the Gortler mode has a significant effect on the 
subcritical breakdown of the TS waves. However, for a sufficiently low level 
of background noise, we should expect that a stationary Gortler mode could be 
set up by slowly increasing the Reynolds number. 

Case (b): Here we consider the interaction of TSL with GL. In this 

situation, other modes of instability can occur at lower Reynolds number but 
since this situation is relevant to the corresponding external boundary layer 
problem we believe it to be of some importance. This is because in this case, 
as in Case (c), the TS wave now bifurcates supercritically and the possibility 
of stable mixed mode solutions must now be investigated. Here we concentrate 
on the interaction of such a mode with a GL vortex. Case (c) will be con- 
cerned with the interaction with a GU vortex. 

The parameters 6 ^ ,6 ^>r\ and ti^ are a H negative and satisfy 

< Bj/B 2 < n 1 /n 2 , (B ± > 0, i = 1,2). 

A routine calculation lead to the bifurcation pictures shown in Figures 
(12a,b). We see that the mixed mode always bifurcates from the "pure-mode" 


which is the least unstable on the basis of linear theory. This bifurcation 
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leaves the pure-mode stable so that at sufficiently large p , both pure 

modes are possible stable equilibrium states. However, in the absence of any 
finite amplitude background noise we expect that the pure mode which is the 
most unstable on the basis of linear theory would be set up when the Reynolds 
number is gradually increased. 

Case (c): Here we consider the interaction of TSL and GU. For this 

case ^i^2 ,r| 2 are ne £ at * ve whilst ri^ is positive. The relationship 
between the ratios of the coefficients is: 

6 1 /6 2 < B 1 /e 2 < -n 1 /n 2 , (S i > 0, i - 1,2). 

The bifurcation pictures are shown in Figures (13a,b). If the Gortler mode is 
the most unstable on the basis of linear theory, then there is no secondary 
bifurcation and the TS mode is never stable. When the TS wave bifurcates 
first, then it is initially stable before it suffers a secondary bifurcation 
to a stable mixed mode. The mixed mode then meets the "pure" Gortler mode 
which changes from being unstable to stable. Thus, for both v > 0 and 
v < 0 at sufficiently large values of p the only stable state possible 

is that corresponding to a finite amplitude Gortler vortex. Hence, the 
Gortler mode effectively prevents the finite amplitude growth of the TS wave. 


6. CONCLUSIONS 

In this paper we have considered the interaction of two types of pertur- 
bations in a curved channel flow; these are the travelling nonaxisymmetric 
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wave (TS) and the axisymmetric vortical perturbation referred to as the 
Gortler vortices. 

By using the Stuart-Watson approach, two coupled equations for the ampli- 
tudes of TS and the Gortler perturbations were obtained. Coefficients of 
these equations have been calculated for their interaction, from Reynolds 
number starting at the common critical value Re c for both the perturbations 
up to a large enough value which we think covers all the possible inter- 
actions. We have, however, concentrated our attention on those interactions 
which we think are significant in external flows. 

We have seen in the previous section that for R close to Re c the only 
possible stable ,f pure state" to be Gortler vortices. For a finite range of 
Reynolds numbers, a mixed mode is possible, but in any experimental investi- 
gation of this problem, we expect this range to be too small to be detected. 

However, the threshold amplitude effect associated with a finite amplitude TS 
wave remains intact and indeed is augmented by the curvature. In external 
flows such as a Blasius boundary layer or an attachment line boundary layer, 
this effect, if repeated, would make these flows more sensitive to background 
noise. 

Consideration of the interaction between a TS perturbation corresponding 
to the lower branch of its neutral curve with a Gortler perturbation belonging 

to the lower branch of its neutral curve, shows that a stable finite amplitude 

perturbation of either type can be set up depending on which one is most 
linearly unstable. The value of the radius ratios q determines which of 
the perturbations is most unstable. 

For the interaction of the TS perturbation corresponding to the lower arm 

•• 

of its neutral stability curve with the Gortler perturbation corresponding to 


-29- 


the upper arm of its neutral stability curve, we find that the Gortler vortex 

prevents the occurrence of a finite amplitude TS wave far from the neutral 

curve. When a TS wave is the most linearly unstable of the two perturbations, 

a finite amplitude TS wave develops, the amplitude of which increases as the 

Reynolds number increases further from its value on the neutral curve until a 

• • 

"mixed" mode appears. Here both the TS and Gortler have finite amplitudes. 
As Reynolds number increases further, the mixed mode bifurcates into a stable 
Gortler mode. In the case when the Gortler mode is the most linearly 
unstable, only a finite amplitude Gortler state is possible as the Reynolds 
number increases from a value on the neutral curve. 

For external flows, an asymptotically self-consistent description of non- 
linear TS waves has been given by Smith (1979). Here the disturbance was 
described by "Triple Deck" theory and the streamwise scaling for the TS wave 
corresponds to lower branch TS waves in our problem. Further it was shown 
that lower branch TS waves bifurcate supercritically so we can expect that our 
results for the interaction of TSL waves and Gortler vortices in channel might 
have implications to the external flow problem. Of course, the effect of 
boundary layer growth might negate the validity of us drawing these con- 
clusions, however, we believe that our calculations show what is the likely 
effect of the possible interactions involving TS waves and Gortler vortices. 

In a later publication, we will report the numerical simulations of the 


interaction in curved channel flow between such perturbations 
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APPENDIX A 


Here the nonzero components of the operators of (3.11) are given explic- 
itly in terms of their row-column location, i.e., (1,2) will refer to element 
in first row and second column. 


I*, [3 ,3„] is a 6x6 matrix: 
1 Z 0 


(1,2) = - 


(1-n) * 

v • 


(1,3) = - 


(1,4) = 


(1-n) 2 „ . 3 zz f(r) (1-n) „ 

2 2 d ee R 2 “ r R 9 

R o r R 0 ° 


(1,5)= - 3 q + | f(r)(l-n) + 3 0 

R 0 r r R o 


( 2 , 1 ) = 


R 0 (l-n) 


( 2 , 2 ) = - 


(1-n) 


(2,4) = - a + (1 - n) + (1 - n) 


R 0 r 


(2,5) - (1 - r,) 2 /r 2 - 3,, - 3^ + liEi (1 - n) 3 ( 

r 


(3,0 - 3 z 


(4,5) » 


R 0 O-n) 
r 

(4,6) = - 3 z 

(5,2) = (6,3) = 1. 


Tj[8 t ] is a 6x6 matrix: 


(1,4) 


ml 3 

0 


(2,5) = (3,6) = R Q (1 - 


L„[3 ,8.] is a 6x6 matrix 
L z u 
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(1.4) 

(1.5) 

(2,1) 

(2.4) 

(2.5) 
(3,1) 

(3.6) 
(4,5) 


R o r R o 


1 

n) ^2 a 0 
R o 


R o r r R o 


■r (1 -" )9 e 


^12 2 
-J-U 3 e 

R o r 


*1^ 


2 R 0 R,3 
0 1 z 


R i < 1 - R > 


-T-d V 


T 2 1 [ ^ t ^ is a 6x6 matrix: 


(1.4) = — i- (1 - n )3 f . 

R 0 

(2.5) = R x (l - n)3 t 

(3.6) = R x (l - n)3 t . 
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T„„[9 ] is a 6 x 6 matrix: 
22 T 


(1,4) = - 


(1-n) -x 

"*r * 


(2,5) = (3,6) = R q (1 - n)3 T 


Ni[3 z ,a 0 ] {<[® <l} is a 6 component column vector: 


(1,1) = -(1 - p) [-£— 3 0 u - + ~~ 2 ~ m f w9 z u ~ u9 z w } - 4 " ^ + u 9 0 V H 


Rj(l - n) 




(3,1) = r q (i - n) y 9 e w + u |^r + w9 z w. 


i 
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APPENDIX B 

3 

Here we give the nonlinear terms of 0(e) which form part of the 
coupled Landau equations. Note the presence of cubic terms such as 

o 

|x(t) | X(x), etc. The lower case letters are elements of the vectors given 

in capitals, e.g., G_ - (gj ,g 2 >g 3 >g 4 >g 5 >g 6 )^« Each of these elements are 
functions of £ . 

(i) Ik • (contribution from N^) = 

(^-) [2j 5 a 5 + 2c 5 a*] X(T)|x(x)| 2 a 1 
+ Oral [ c ^ a * + c 5 a* + j 4 a 5 + j $ a 4 ] X(x ) |x(x ) | 2 a 2 
+ [-ik c^ a $ + 2ik a & c $ + j & ika 5 ] X(t ) | x(r ) | 2 a 2 

+ [ c 4 a 2 + a 4 c 2 + j 4 a 2 + a 4 j 2^ X(t > l x < T > I \ 

+ [~ik c^ a* + 2ik a* c & + j^ikag] X(x) |x(x) | 2 aj 
R 0 

+ [c 4 a* + j 4 a 4 ] X(x)|x(x)| 2 a^ 

R o r 

- [~ik c 6 a* + 2ik a* c 4 + j^kaj X(x ) | X(x ) | 2 Ij 
R 0 

+ [(c 4 a* + a*c 3 ) + (j 4 a 3 + a^)] X(x) |x(x) | 2 a 3 


+ [-ik c 6 a & + 2 ik a & + jgi^] X(x)|X(x)| Z a 3 


(1-n ) 


4^- [ _im h 4 b 4 + im b 4 + g 5 im b 4 “ ira b 5 S 4 ] X(t)|y(t)| 


+ [ b 6 ik h 5 + b & ik g $ + ik a^] X(x ) | Y(x ) | 2 a 2 


Rq -- 1 ' ^ ~ [~im b 5 t > 5 + b $ im h $ + g 5 imb 5 -im b $ g 5 ] X(x)|Y(x)| 


* * 


+ U_Hl [h 4 b* + b*h 5 + g 4 b 5 + b 4 g 5 + a 4 k 5 + a^] X(x ) | Y(x ) | 


* * 


+ [h 4 b 2 + b"h 2 + g 4 b 2 + b 4 g 2 + a 2 k 4 + a 4 k 2 ] X(x)|Y(x)| 2 a 2 


R Q (l-n) 


* * 


- [-im h $ b & + b $ im h & + g $ im b g -im b $ g & ] X(x)|Y(x)| 


+ lilUl [2(h 5 b 3 + g 5 b 5 + a 5 k 5 )] X(x)|Y(x)| z l 1 


+ [b 6 ik h 6 + b 6 ik g fi + k & ik a & ]x(x ) |Y(x ) | 2 a 3 


+ [h.b* + b*h + g b + b g + a,k. + a.k ] X(x ) | Y(x ) | 2 a, 


+ [-im h 4 b* + b* im h 5 ] X(x ) | Y(x ) | 2 a 1 


- [b*ik h 4 + b 6 ik g 4 + ik a 4 k 6 ] X(x)|Y(x)| Z ^ 


+ ~ 2 ~ [ b 4 b 4 + g 4 b 4 + a 4 k 4 ] X(x)|Y(x )| 2 
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+ iy [b* ik h 6 + b 4 ik g 6 + ik a 6 k 4 ] X(x)|Y(T)| 2 a 1 
R 0 

+ t 1 " 1 g 4 b 5 ' lm b 4 g 5 ] X ( T >I Y(T) I% 

(ii) J5 T • (contribution from N^) = 

+ llril [2 d 5 b* + 2k 5 b 5 ] Y(T)|Y(x)| 2 b 1 

_ [-im d 5 b* + 2im b* d 4 + k 5 im b 4 ] Y(x ) | Y(x ) | 2 bj 

+ [d 4 b* + b* d 2 + k 4 b 2 + b 4 k 2 ] Y(x)|Y(x)| 2 b 2 
R n ( 1 ^ ^ 

+ — f-im d c b_ + 2im b c d c + k- im bJ Y(t) Y(t) b 0 

r L 55 5 5 5 5 J l 

* [ d 4 b 5 + d 5 b 4 + V 5 + b 4 k 5 l «t)|Y(t)|*b 2 

+ - — [-im d^ bg + 2im b^ d^ + k^im b^] Y(x)|Y(x)| b^ 

+ [d 4 b* + b*d 3 + k 4 b 3 + b 4 k 3 ] Y(x)|Y(x)| 2 b 3 

+ 2r l 2 [-im d 4 b* + 2im b* d $ + k 5 imb 5 ] X(x ) | Y(x ) | 2 b 1 

+ 2(1 rUl [d 4 b* + k 4 b 4 ] Y(x)|Y(x)| 2 b 1 

R o r 

+ — [^ 5 a* + g*a 5 + b 5 j 5 ] Y(x)|x(x)| 2 b 1 
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“ [a 5 im h 4 + a 5 im g 4 + b 4 imj 5 ] Y(t ) |x(t ) | 2 % l 

+ [\ a * 2 + a 4 h 2 + zl a 2 + a 4 g 2 + b 2 j 4 + b 4 j 2^ Y(x > • X(t ) I \ 

+ -^1 R q [ a* imh 5 + a 5 im g* + imb^] Y(x ) |x(x ) | 2 b 2 

+ [h 4 a 5 + a 4 h 5 + g 4 a 5 + a 4 g* + b 4 j 5 + b^J Y(x) |x(x) | 2 b 2 

+ [-ik h & a 5 + a 6 ik h 5 + g*ika $ -ik a & g*] Y(x ) |x(x ) | 2 t> 2 

(l-n)R Q * * 9 

[a 5 im h 6 + a 5 im g & + imj^] Y(x)|x(x)| Z b 3 

+ [h 4 aj + 8484 + b 4 j 4 ] Y(x)|X(x)| 2 b 1 

R o r 

+ [ a 4 im h 5 + a 4 im g* + inrt^jj Y(x ) | X(x ) | 2 bj 

+ ~Y [“ lk h 4 a 6 + a 4 ikh 6 + g 4 ika 6 " ik a 4 Y(x) I x (t) 1 2 bj 

R 0 

- [" ik h 6 a 4 + a 6 ik h 4 + g 6 ika 4 “ ik a 6 g D Y(x > l x < T ) | 2 bj 
R 0 

+ [h 4 a* + a*h 3 + g*a 3 + a 4 g* + b 3 j 4 + b 4 j 3 ] Y(x) |x(x) | 2 b 3 . 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


FIGURE CAPTIONS 


Curved channel with flow in the azimuthal direction. The walls 

are parallel to the z-axis. Radius ratio n = R. /R and 

1 o 

channel width d = R Q - R^. 

Critical Reynolds number versus (1 - n) for Gortler and TS 

perturbations. At the cross-over point of these curves. 

Re - 8 x 5772.2 and (1 - r\ ) = 2.179 x 10~^. This figure 
c c 

is a schematic adaptation of Figure 1 in Gibson & Cook (1974). 

Neutral stability curve for the Gortler perturbation. Reynolds 
number versus axial wave number k. 

Neutral stability curve for the Tollmien-Schlichting perturba- 
tion. Reynolds number versus azimuthal wavenumber m. 

*• 

The two neutral stability curves (a) Gortler neutral stability 
curve (b) TS neutral stability curve for n = n c » For this 

case, the critical Reynolds number is identical for the two per- 
turbations. GL and GU refer to the lower and upper arms of the 
Gortler stability curve. TSL and TSU refer to the lower and 
upper arms of the Tollmein-Schlichting stability curve. 

Reynolds number versus growth rate 0^ of the Gortler pertur- 
bation. 

Reynolds number versus growth rate 0 ^ of the TS perturba- 

tion. Note the negative 8^ for P art of TSU. 

Reynolds number versus for perturbations corresponding to 

the two arms GL and GU of the Gortler stability curve interacting 
with TSL and TSU: 
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(a) GL interacting with TSL and TSU. 

(b) GU interacting with TSL and TSU. 


Figure 9 Reynolds number versus 


for perturbations corresponding to 


the two arms GL and GU of the Gortler stability curve interacting 
with TSL and TSU: 

(a) GL interacting with TSL and TSU. 

(b) GU interacting with TSL and TSU. 

Figure 10 Reynolds number versus r^* 

• • 

Figure 11 Bifurcation diagrams for the interaction of the Gortler and TS 
perturbations for Reynolds number at or close to Re c where 


n j /n 2 > > ^l^ 


n i > 0, B i > 0, < 0, i = 1,2. 


(a) (v > 0). The Gortler mode is the most linearly unstable 
mode in this case. The TS mode is subcritically unstable. 

(b) (v < 0). The values of y at P and Q are 


respectively 


- v V 6 i 

n Bo 

11 


-'V S l 

5 , 6? 

(I 


Figure 12 Bifurcations diagrams for the interaction of GL and TSL where 


n 1 / n 2 > e 1 / 6 2 > $ l /& 2 
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ri i < 0, 8 ± > 0, 6 < 0, i = 1,2. 


(a) (v < 0) . The TS mode is the most linearly unstable mode in 
this case. The value of y at Q is 




- v V b i 

6, e. 


(b) (v >0). The Gortler mode is the most linearly unstable 
mode in this case. The value of y at P is 


y P = 


-V B i 

n, e 7 

(1 - J. JL) 
V "2 V 


Figure 13 Bifurcation diagrams for the interaction of GU and TSL where 


and 


- h 1 /h 2 > > 


> 0, n 2 < o, > 0, 6 i < 0, i - 1,2. 


(a) (v < 0). The TS is the most linearly unstable mode in this 
case. The values of y at P and Q are respectively: 


y P = 


and 


- v V 8 i 

n 1 


B q‘ 


- v V e i 


(b) (v > 0) • The Gortler mode Is the most linearly unstable 
mode in this case# 
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